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Abstract. The response function of a network of springs and masses, an 
elastodynamic network, is the matrix valued function W(u), depending on the 
frequency ui, mapping the displacements of some accessible or terminal nodes to 
the net forces at the terminals. We give necessary and sufficient conditions for a 
given function W(oj) to be the response function of an elastodynamic network, 
assuming there is no damping. In particular we construct an elastodynamic 
network that can mimic a suitable response in the frequency or time domain. 
Our characterization is valid for networks in three dimensions and also for 
planar networks, which arc networks where all the elements, displacements and 
forces are in a plane. The network we design can fit within an arbitrarily small 
neighborhood of the convex hull of the terminal nodes, provided the springs and 
masses occupy an arbitrarily small volume. Additionally, we prove stability of 
the network response to small changes in the spring constants and/or addition 
of springs with small spring constants. 



1. Introduction 

Is it possible to design an elastic material that has a prescribed response? This 
question is answered by Camar-Eddine and Seppecher Q for linear elastic materials 
in three dimensions, assuming the macroscopic response is governed by a single 
displacement field. Their approach consists of three steps. First it is shown how 
to design a continuum material that behaves like an elastic network (a network 
composed of springs). Then the response of elastic networks is characterized, i.e. it 
is shown how to construct an elastic network with a suitable response. A limiting 
argument is then used to answer the question for the continuum. As a first step 
towards solving the characterization problem when the response depends on time, 
we show how to design an elastodynamic network (a network of springs and masses), 
that can mimic a prescribed response as a function of time (or frequency) . Moreover 
if the springs and masses occupy an arbitrarily small volume, the network can 
be designed to fit within an arbitrarily small neighborhood of the convex hull of 
the terminal nodes, which is a requirement for an argument similar to that of 
Camar-Eddine and Seppecher 0]. An earlier characterization of elastodynamic 
networks is that of Milton and Seppecher . However the network elements used 
in the construction || are frequency dependent, so the constructed network can 
only mimic the response function at a single fixed frequency. 

In a different context, the approach of Camar-Eddine and Seppecher was ap- 
plied earlier by the same authors Q to characterize all possible responses for the 
conductivity equation, assuming the macroscopic response is governed by a single 
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voltage field. The problem of finding a network with a given response is often called 
"network synthesis" , and the earliest example is Kirchhoff's Y — A theorem, which 
characterizes the response of any resistor network in three dimensions. Another 
characterization for resistor networks is that of Curtis, Ingerman and Morrow [B[ 
who consider planar networks that can be embedded inside a disk and where all 
terminals lie on its boundary. For electrodynamic networks (with resistances, ca- 
pacitors and inductances), we are only aware of results dealing with the frequency 
response or impedance of a circuit with two terminals (see Foster 0, 0] and Bott 
and Duffin Milton and Seppecher Q give a construction for n— terminal elasto- 
dynamic, electrodynamic and acoustic networks which is valid at a single frequency. 
The electro mag netic analog of elastodynamic networks is considered by the same 
authors 1(J mi 

In |j2]we give the properties of the response function of elastic and elastodynamic 
networks. The construction of a network that matches a response function with 
the properties in |J5]is given for the static case in 33 Note that the characteriza- 
tion of elastic networks by Camar-Eddine and Seppecher Q is part of a limiting 
argument on energy functionals, so only non-degenerate three dimensional clastic 
networks are needed. The degenerate case corresponds to planar elastic networks 
(the network, forces and displacements lie on a plane) and is a set of measure zero 
which leaves the energy functionals considered in Q unaffected. We complete the 
characterization in Q to include planar elastic networks. Then in Sj4] we completely 
characterize the response of elastodynamic networks (planar or in three dimensions) 
for all frequencies and assuming there is no dissipation (damping) in the network. 
We include in the appendices two technical results. Appendix [A] shows that the 
network response is stable with respect to small changes in the spring constants 
and the addition (but not deletion) of springs. Appendix [B] uses stability to give 
a systematic method of modifying an elastic network to eliminate floppy modes 
without changing significantly the response. Floppy modes correspond to nodes 
that can move with zero forces and they are discussed in more length in 



1.1. Preliminaries. Consider a network composed of springs and masses, and 
assume we only have access to n "terminal" or "boundary" nodes Xi, . . . , x„ £ M. d , 
where the dimension d is either 2 or 3. The network is said to be planar if d = 2 
and the springs do not cross. The static response matrix or displacement-to-forccs 
map is the nd x nd matrix W so that 

f = Wu, 

where u = (u[ , . . . , uJ l ) T is the vector of displacements of the terminal nodes Xj 
and f = (ff, . . . , ) T is the vector of net forces acting on node x» at equilibrium. 

In the dynamic case the displacements u(t) and f (t) depend on time t. Let u(oj) 
be the Fourier transform of u(i), 



/oo 
u(t)t 
-oo 



-iujt 



dt. 



and similarly for f(w), where w is the frequency. Then if lu is not a resonance 
frequency of the network (a precise definition of resonance is given later in fl2.2.2D , 
the response matrix of the network is the possibly complex nd x nd matrix valued 
function W(w) such that 

f (w) = W(w)u(w). 
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For convenience we have chosen to work in the frequency domain. However when 
u(t) = for t < 0, our results can be reformulated for the transfer function of the 
network since £[u(t)](s) = u(— is), where C denotes the Laplace transform, i.e. 

/•OO 

£[u(t)](s) = / u(t) e - st dt. 
Jo 

In this case the transfer function of the network is W(— is). As we work only in 
the frequency domain, we drop the hats in the Fourier transform notation for the 
sake of clarity (i.e. u(cu) = vL(u>) etc. . .). Also as there is no dissipation, it suffices 
to assume that u(w) and f (u>) are real to determine the real valued function W(w). 



2. The response function of an elastodynamic network 



In this section we establish the properties that the response of an elastodynamic 
network satisfies. We start with the response of networks (static or dynamic) where 
all the nodes are terminals ( £|2.ip and then study the case where interior nodes are 
present ( §2.2p . We also include some transformations in §2.31 that do not affect the 
response function. 

2.1. Response function for networks without interior nodes. Consider the 
simple network consisting of two nodes xi and X2 with masses m\ and m 2 , linked 
with a spring with spring constant fci,2- Let be the force exerted by the spring 
on node Xj, i = 1, 2. By Hooke's law 

(x 2 - xi)(x 2 - Xx) T 

a 2 = -k Xft -2 (u 2 - ui) = -ai. 

11*2 - xi || 

The laws of motion can be written in matrix form as — cj 2 Mu = f — Ku, where 



K = fci a 



ni j2 n^ 2 



-ni i2 n^ 2 
n i,2nf 2 



M = diag (mie, m 2 e), 



ni, 2 



x 2 - Xl 



x 2 



and the vector e = (1, . 
single spring is given by 

(1) 



■xi || 



for d = 2,3. Thus the response function of a 



W(w) = K uj 2 M. 



When all nodes are terminal nodes (i.e. there are no interior nodes) the response 
function can also be written in the form |T]), but now K is the stiffness matrix of 
the network and the mass matrix M = diag (m\e, . . . , m n e), where m,; > is the 
mass of the i— th node and n is the number of nodes. The stiffness matrix of the 
network is the sum of the stiffness matrices associated with the individual springs, 



vndxnd 



springs i,j 



where the summation is over all pairs of nodes x, ,Xj connected by a spring. The 
2d x 2d matrix H£(i j) is the stiffness matrix of the spring between nodes Xj and 
Xj. We also used nd x d matrices Ej = [e(j_xw+i, ■ • ■ > which are introduced 
so that the components of 'K.(i.j) enter the appropriate blocks of K (e p is the p— th 
canonical vector in M. nd ). We only consider non-negative spring stiffnesses k 



hi ■ 
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Stiffnesses with a non-zero imaginary part model damping or dissipation of energy 
in the network and are left for future studies. 

2.2. Response function for networks with interior nodes. 

2.2.1. The static case. The response matrix W can be obtained from the response 
matrix A of the network where all nodes are considered as terminal nodes. The 
partitioning of the nodes into interior nodes / and terminal (boundary) nodes B 
induces the following partitioning of A, 

(2) A = \ A A BB K K m . 

v ' [A IB A u 

Instead of dealing directly with the response matrix A, it is convenient to introduce 
the quadratic form 

9a(u) = u T Au, 

which represents twice the total elastic energy stored in the network. In the simple 
case of a single spring between nodes Xi and x 2 with spring constant k, the quadratic 
form is 

xi - x 2 



S( Xl .x 2 )(ui,u 2 ) = k (ui -u 2 ) 

V ll X l _ X 2| 

We omit the spring constant indices for clarity. When there are more springs qx is 
the sum of the quadratic forms for all springs, thus <?a(u) > 0. 

For general static networks the response matrix is defined indirectly by its qua- 
dratic form q-w: 

(3) (?w(u B ) = inf g A (us,u/). 

By the partitioning (f2]) we may rewrite 

(?a(ub, u/) = UsA BB u B + 2115AB/U/ + uj A//U/. 

The first order optimality conditions for the minimization (J3j> are actually the bal- 
ance of forces at the interior nodes: 

= V U/ q A (u B , Uj) = 2 A//U/ + 2 A/ B ub • 

The following lemma shows that for any Ug it is possible to balance forces at the 
interior nodes and it implies the minimization has at least a minimizcr (since 
</a is bounded below). Another way of seeing this lemma is that if there are any 
"floppy" modes within the interior nodes (i.e. modes that generate displacements 
with zero forces) then those modes are not coupled to the terminals. 

Lemma 1. Given the partitioning (f2]) of the response matrix A where all nodes 
are considered as terminal nodes, we have IZ(Aib) C TZ(Ah). Here 7?.(B) denotes 
the range of a matrix B . 

Proof. By reciprocity A T = A. thus it is equivalent to prove AT(Abi) 3 AT(An), 
where A/"(B) denotes the nullspace of a matrix B. Let u/ be a displacement such 
that A//U/ = (i.e. a "floppy" mode). Then 





= uj A//U/ = [0 uj] A 



springs i, j£l springs i£l,j£B 
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Therefore for all nodes x, £ I and Xj £ B that are linked by a spring we must have 
u; • (x, — Xj) = 0, which means precisely that Ab/u/ = 0. □ 

Remark 1. We show later in Appendix \B\ that floppy modes can be eliminated 
from a network by adding springs with small spring constants. The response of the 
new network can be made arbitrarily close to that of the original one, provided the 
new springs have sufficiently small stiffness. Examples of floppy modes are given 
in Figured 

By eliminating the interior nodes, the static response matrix can thus be written 
in Schur complement form: 

(4) W = A B b-A bj aJ j Ajb j 

where f stands for the Moore-Penrose pseudo-inverse, which is is simply the inverse 
if there are no floppy modes. 

We denote by uA v the cross product of the vectors u, v £ M. d . For d = 2 we have 
u A v = det[u, v] and for d = 3, u A v = (U2V3 — U3W2, U3V1 — U1V3, u\V2 — U2V\) T . 
Before reviewing some properties of the static response matrix we need the following 
definition. 

Definition 1. A balanced system of forces fj, i = 1, ...,n supported at nodes Xj, 
i = 1, . . . ,n in M. d (d = 2, 3) satisfies: 



(a) /J fj = (balance of forces) 

i=l 
n 

(b) ^""^ Xj A f ; = (balance of torques) 



i=l 

Lemma 2. The static response matrix satisfies the following properties. 

(a) W £ W ndxnd . 

(b) W = W T (reciprocity) 

(c) W is positive semidefinite (energy is not produced by the network) 

(d) Every column f = (f^, . . . , f^f) T of W is a balanced system of forces when 
supported at the nodes Xj in R d . 

Proof. Properties (a), (b) and (d) follow from the construction of the response 
matrix. We now prove Property (c). Let W be the response matrix of a network if 
all the nodes are considered as terminal nodes. Then for all displacements u £ M. nd 
we have 

q(u) = u T Wu = ^2 ^^(u^Uj) > 0, 

springs i,j 

where S( x . , Xj )(u, v) is the quadratic form associated with the spring between nodes 
Xj and Xj. Thus (c) holds for networks where all the nodes are terminals. Using 
([3]) wc sec that (c) holds for general networks as well. □ 

2.2.2. The dynamic case. The response function in the dynamic case can be ob- 
tained in a similar way as in the static case. First if all the nodes are terminal 
nodes, the response function A(w) of the network is given by |T]). The partitioning 
of A induced by the partitioning of the nodes into boundary B and interior / nodes 
is, 



AH 



Kbb Kb/ 



Mbb 



(i 
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As in the static case we can introduce the quadratic form 

(5) g A (u B! U7;w) = u^(K B b - uj 2 M bb )u b + 2u|K B/ u/ + uj (K n - w 2 M//)u/. 

Remark 2. Unlike in the static case the quadratic form qA(u B , U/; could be 
unbounded from below for u B fixed. This happens for example if there is a u/ so 
that K//U/ 7^ and M//U/ ^ 0. Then for lj large enough the matrix in the last 
term of qx becomes indefinite. Thus we cannot define the response function at the 
terminals through a minimization principle similar to (|3|). 

The dynamic response function at the terminals is the displacement-to-forces 
map at the critical point V U/ (?a(ub, u/; uj) = 0, if such critical point exists. Because 
<7A may be unbounded below, this critical point could be a saddle point for the 
quadratic <ja with u B fixed. The frequencies u> for which there is no critical point 
(i.e. there is some u B so that V U/ <7a(u.b, u/; uj) ^ for all u/) are important 
physically and in our derivation and correspond to the resonance frequencies of the 
network. 

To give an expression for the dynamic response function we partition the interior 
nodes into nodes J with positive mass and massless nodes L, so that I = J U L. 
Therefore Mjj is positive definite but M^l = 0. 

Lemma 3. The response function at the terminals is 

(6) W(w) - K BB - uj 2 M bb - K bj (Kjj - w'MjjJ^Kjb, 

provided that uj 2 is not an eigenvalue o/Mjj^KjjMy]' 2 . Here we have used the 
submatrices of the matrix 

it\ xr K BB K B j K BB K BJ K Bi „t [i- v 1 

K JB Kjj [>^JB K-JjJ L Kji J 

Proof. The matrix K is the response matrix for the network with terminals B U J 
and interior nodes L and can be obtained from ^ . Since the nodes L are massless 
the dynamic response at the nodes BU J is K — cj 2 diag (M. BB , Mjj). Since Mjj is 
non-singular, the matrix Kjj — oj 2 M,/j is singular if and only if uj 2 is an eigenvalue 
of M}y 2 KjjMjj /2 . Thus when uj 2 is not an eigenvalue of Mjj /2 KjjMjj /2 , 
we can equilibrate forces at the nodes J and get the expression for the response 
function. □ 

A corollary of Lemma [3] is that if lj is a resonance frequency of the network then 
to 2 must be an eigenvalue of the matrix Mjj^KjjMjj^ 2 . The expression for the 
response function in Lemma [3] leads to the following properties. 

Lemma 4. The response function W(o>) of any network of springs and masses 
with n terminals is of the form 

(8) W(uj) = A - uj 2 M + V — =• e w idxnd , 

i uj uif 

where the matrix M = diag (mie, . . . , m n e) is real diagonal with the masses of the 
boundary nodes in the diagonal, the vector e = [1, . . . , 1] T G R d , the matrices 
are real symmetric positive semidefinite, and the static response 

W(0)=A-5>r 2 C«, 
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is real symmetric positive semidefinite and balanced (i.e. it satisfies the conditions 
(a)-(d) of Lemma\^). The resonant frequencies are distinct, finite and satisfy 
co 2 > 0. 

Proof. Let I = J U L be a partition of the interior nodes into massless nodes L 
and nodes with positive mass J, and let K be defined as in ([7]). By Lemma [3] the 
response function at the terminals can be rewritten as 

W(w) = K BB - lu 2 M bb - K BJ M7j /2 (C - u^M^Kje, 

where C = M~j /2 KjjMjj /2 and ui 2 is not an eigenvalue of C. The matrix C 
is symmetric positive semidefinite because K is a symmetric positive semidefinite 
matrix (LemmaEJ). Let {oj 2 }f =l be the (nonnegative) eigenvalues of C and {cjjjLi 
be a corresponding orthonormal basis of eigenvectors of C, where N = \J\. When 
u ^ and uj 2 ^Wj, the response function W(w) becomes 

N ~~T 

(9) WH = K BB - lu 2 M bb + V J " 3 2 , 

— 1/2 — 

with Cj = K B jMjj Cj, j = 1,...,N. Let r = rank(C) = rank(Kjj) and 
assume the eigenvalues are ordered such that u> 2 > for j = l,...,r. Clearly 

KjjMj] /2 z = if and only if Cz = 0. Thus Mjj /2 Cj G Af(Kjj), for j = 
r + 1, . . . , N. By Lemma [TJ we have TV'(Kjj) C JV(Kj B ) which means that 

Zj = for j =r + 1, ...,N. 

In other words, only the first r terms of the sum in @ are nonzero. We obtain the 
form ([5]) of the response function from ^ by setting A = J£ BB and M = M. BB . 
The matrices are the sum of the matrices Cjcj that correspond to the same 
resonance oj 2 , thus the must be real symmetric positive semidefinite. 

We now show that W(0) = J£ BB — K B jK\jKj B , i.e. at = the dynamic 
response function is the static response of the network. Then the properties of 
W(0) follow from Lemma [5] First note that from 

W(0) = K BB - K BJ Mjj /2 &Mjj /2 Kj B 

where we used that 

It is sufficient to show that u,/ = — Mjj^C^M. J 1 / 2 'Kj B u B equilibrates the forces 
at the interior nodes for any terminal displacements u B . Indeed we have 

(10) CM^uj = (Mjj^KjjMjj^MYjUj = -M7] /2 K JB u B , 
since Lemma Q] and Mjj invertiblc imply 

11(C) = Tl(M-j /2 Kjj) D K(M-y 2 K JB ). 

The balance of forces at the nodes J (i.e. KjjUj = — K.j B u B ) comes from multi- 
plying (ED by M 1 / 2 on the left. □ 
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Figure 1. Truss structure replacing a spring between nodes Xj 
and Xj without changing the response function. This structure 
can be used to avoid a line for (a) networks in R 3 or (b) planar 
networks. 



2.3. Network transformations not affecting the response function. We 

need a few elementary transformations that allow more flexibility with the place- 
ment of springs in a network. We assume throughout this text that the springs 
occupy an arbitrarily small volume and that the nodes are points which may or 
may not have a mass attached to them. 

2.3.1. Avoiding a line. It is possible to transform a spring in order to avoid a line 
or a crossing (for networks in R 3 ). The construction for networks R 3 is given in 
({J Example 3.2] and consists of replacing the spring by a simple truss as shown 
in Figure [lja). A similar construction can be done for planar networks, see Fig- 
ure [ljb). If a network in R 3 has springs crossing, these can be eliminated via 
this transformation since the position of the additional interior nodes u, v, w is not 
fixed. Moreover, the additional nodes in the truss structures can be chosen to avoid 
a finite number of points. 

2.3.2. Virtual crossings. A network with all springs in R 2 is not necessarily planar 
because its springs may cross. However [9j, Example 3.15] shows how to replace 
such a crossing by a planar network with exactly the same response function. This 
transformation involves adding a node at the crossing point of the springs and 
carefully choosing the spring constants. To avoid a finite number of points one can 
first replace one of the springs by a simple truss as in Figure [TJb) and use virtual 
crossings to transform the network into a planar network. 

2.4. The superposition principle. A fundamental tool for our construction of 
a network reproducing the response function is the following result, which is valid 
for both planar and R 3 networks. 

Lemma 5. Let Wi and W2 be the response matrices of two networks (planar or 
in R 3 , static or dynamic) sharing the same terminals but with no interior nodes in 
common. Then the response function of both networks together is Wi + Wj . 

Proof. The result follows from the reasoning in [^, Remark 3.9] and the frequency 
independent transformations in ^2.31 For the planar case any crossing can be elim- 
inated using [^, Example 3.15]. Note that the transformations in i j2. 31 allow one to 
avoid a finite number of locations (except the terminals). □ 
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3. Characterization of the static response 

Building upon the seminal work of Camar-Eddinc and Scppecher we give 
necessary and sufficient conditions for a function to be the response function for 
either a planar network or a network in K 3 . The sufficiency is proved constructively 
and relies on the existence of networks that have rank one response matrices, as is 
described in detail in the remaining part of |J3] 

Recall that the e-neighborhood C e of a set C is the set, 

(11) C £ = {x e 1 2 dist(x, C) < e}. 

If the set C is convex then the set C € is also convex because of the convexity of the 
function dist(x, C) for convex C (see e.g. @, §3.2.5]). 

Theorem 1. For any choice of terminal node positions, any matrix W satisfying 
the properties in Lemma\^ is the response matrix of a purely elastic network which 
is either planar or in R 3 . Moreover, any internal nodes in the construction can be 
chosen within an e— neighborhood of the convex hull of the terminals, and avoiding 
a finite number of positions. 

Proof. By properties (a)-(c) in Lemma [2] the matrix W can be written as a sum 
of rank one matrices W.;: 

n 

4=1 

where Wj = AjW^w^ and (A^w^) is an eigenpair of W, Aj > 0, for i = 1, . . . ,n. 
Each satisfies properties (a)-(d) in Lemma [H Properties (a)-(c) are easy to 
check for W and (d) follows by linearity, since it holds for each Wj. Owing to 
Theorem [3] (Theorem [5] in the planar case) , it is possible to construct a network 
with matrix response equal to W.;. By the superposition principle we obtain a 
general network with response W. If the desired network is planar, then every 
crossing between springs can be transformed to a planar network through a truss- 
like structure Examples 3.2, 3.15]. As discussed in 32. 3( such transformations 
can be chosen to avoid a finite number of points in R d . □ 

Remark 3. // W is the static response matrix of a network and a is a positive 
constant then clearly aW is the response matrix of the same spring network, but 
where all the spring constants are multiplied by a. 

3.1. Planar networks with rank one static response matrices. The main 
result in this section is Theorem [2] which is the statement of TheoremQ]for rank one 
response matrices, i.e. it shows that for any rank one response matrix satisfying 
Lemma [2] it is possible to find a planar network that realizes it. We first prove 
Theorem [5] for three terminal networks in 33.1.11 then for four and more terminals 
in 33X21 

3.1.1. Three terminal rank one static planar networks. We show how to construct 
a three terminal planar elastic network realizing any valid rank one response ma- 
trix. If fo,fi,f2 G R 2 is the balanced system of forces at the nodes xo,xi,X2, the 
construction depends on rank [fi, ¥2, Xi — Xo,X2 — xo]. More precisely Lemma [7] 
corresponds to the case when this rank is two and Lemma [8] when this rank is 
one. Since we are in R 2 these are the only non-trivial cases available, which shows 
Theorem [2] for planar three terminal networks. 
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Figure 2. A planar network with rank one response. Any result- 
ing force at terminals is proportional to (f^f, f^, fJ) T . 

Remark 4. (Pierre Seppecher, private communication) The easiest way to con- 
struct a three terminal rank one network is to add a node at the intersection of the 
force lines ( three forces that are balanced meet at a single point in 2D, this can be 
shown by writing the torque balance equation for the intersection point). The only 
problem with this construction is that the extra node can end up far away if the 
force lines are almost parallel. 

We start with the following intermediate result. A similar result is shown in 
three dimensions by Camar-Eddine and Seppecher 0, Lemma 5]. 

Lemma 6. Let fo, fi, fb be a set of balanced forces at the nodes xo, Xi, X2 in M 2 . 
Then if rank [fi, f2, xi — xo,X2 — xo] = 2, there is an e > such that the points 
Xo, x^ = xi + efi and x 2 = X2 + ef^ are not collinear. Moreover e can be chosen 
arbitrarily small and so that x' x and x 2 do not coincide with a finite number of 
points. 

Proof. If it were true that for all e > the three points xo, x'i and x 2 are collinear, 
then the second degree polynomial in e, p(e) = det(x' 1 — xo,x 2 — xo) is identically 
zero. That the constant coefficient of p(e) vanishes means that det(xi — Xo,X2 — 
xo) = 0, or that the points Xo, Xi and X2 are collinear. Since the lemma for 
xo = xi = X2 is trivial to prove, we may assume without loss of generality that 
there is some ael such that 

(12) x 2 - x = a(xi - x ), 

swapping the indices 1 and 2 if necessary. Since the coefficient in e of p(e) vanishes 
we get det(fi, x 2 — x )+dct(xi — x , fb) = or equivalently det(xi — x , — afi+f2) = 
0. Now the torque balance implies that 

(13) (x 1 -x )A(fi+af 2 ) = 

Putting both (|12|) and (fl"3|) in matrix form, there are some real (3 and 7 such that 



I at 




"fi" 




'/3(xi 


- x ) 


-al I 




A 




7(xi 


- x )_ 



The determinant of the matrix above is (a 2 + 1) 2 7^ 0, thus rank [f lt f2,Xi — x , x 2 — 
xo] = 1 which contradicts the hypothesis of the lemma. Finally since p(e) is not 
identically zero, one can choose an arbitrarily small e that avoids a finite number 
of points. □ 

Lemma 7. Let fo, fi, f2 be a set of balanced forces at the nodes xo, xi, X2 in 
M 2 . // rank [fi , f2 , xi — xq,X2 — xo] = 2 then there exists a purely elastic planar 
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network with force response proportional to (f^f, f^, f J) T . The internal nodes of 
such a network can be chosen to avoid a finite number of points and within an 
e— neighborhood of the convex hull of the terminals. 

Proof. First observe the hypothesis that rank [fi , f 2 , xi — xo,X2 — xo] = 2 implies 
the existence of at least one permutation a of {0, 1, 2} such that 

(14) f<r(l) A(x CT(1) - Xo . (0) ) ^0 
Indeed, if (fl4|) is false for all permutations a, i.e., 

fi A (xi - x ) = = f 2 A (x 2 - x ), 

h A (x 2 - xi ) = = f A (x - xi), 

fo A (x - x 2 ) = = fi A (xi - x 2 ) 

we have that rank [fi , f 2 , xi — Xo,x 2 — xo] = 1 and this contradicts our initial 
hypothesis. So, without loss of generality we can assume that, 

(15) fi A ( Xl - x ) ^ 0. 

Next, let x' ; = x.; + efi, for i = 1, 2. By Lcmma[l)J if rank [fx, f 2 , xi — xo, x 2 — xo] = 
2, there is e > such that the points xo, x^ and x' 2 are not collinear and do not 
coincide with a finite number of points. Consider the network in Figure [5] (the 
spring constants are irrelevant for this proof). Let A be the response matrix for 
the network including both terminal (B = {0,1,2}) and interior (I = {l',2'}) 
nodes. The associated quadratic form is 

2 2 

(16) ?a(ub,u/) =^s (x0)X / ) (uo,u^+^s (x . i ^ ) (u i ,u^) + s (x / )Xi) (ui,u^). 

1=1 1=1 
It suffices to show that the quadratic form ^m(us) at the terminal nodes has 
codimension one (since M is positive semidefinite, we do have kerM = kergjyi)- 
Actually = (u^, uf , u^) T G ker^M if and only if there exists uj = (u'^ ,u 2 T ) T 
such that all terms in the sum (] 16|) vanish or equivalently, 

(17) (u'i - u ) T (x' i - x ) = 0, for i = 1, 2, 

(18) {u'i - Utft, = 0, for i = 1,2, 

(19) (u 2 -u' 1 f(x^-x' 1 )=0. 
Property (fT7|) is equivalent to 

u[ - u = oRi(x' - x ) 
u' 2 - u = 6R ± (x' 2 - x ) 

for some a, b reals and where 

" 1" 

-1 ' 



Ri 



Property (|19[) implies that b = a, i.e. the only infinitesimal deformation which 
does not change the side lengths of the triangle in Figure [2] is an infinitesimal rigid 
motion (translation plus rotation). Since f; Ax- = fi Axj and the forces and torques 
are balanced we conclude from (TT5)) that 

222 
^ f, • - u ) = a ^2 f, A (x- - x ) = a f t A Xj = 0. 

i=l i=l i=0 
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Thus ker^M C span {(f T , ff, f 2 T ) T } ± . 

To prove the other side of the inclusion we start with u,, i = 0, 1, 2 satisfying 

2 

(20) ^fi- (Ui-uo) =0 

i=l 

and we seek uj, i = 1,2 such that (fT7|) . (fl8|) and (p~9|) hold. For any real a, the 
choice 

u' : = u + aRi(x'j — x ) and u 2 = u + aRj_(x 2 — x ), 

satisfies (fl7|) and (Tl9|l . Next, note that from the definition of the points x' 1; x 2 we 
have, 

(x- - x ) A fj = (xj - x ) A f i7 for i = 1, 2. 
Using the latter, the balance of forces and (fl"5|) . property (fT8|) follows by taking 

a = fi ■ (ui - u ) = f 2 ■ (u 2 - u ) 
fi A (xi - x ) f 2 A (x 2 - x ) ' 

This proves that kergM = span {(fg , ff, fj) 7 "} -1 , so can be written for some 
c > as, 



□ 



Lemma 8. Let fo, fi, f2 6e a sei of balanced forces at the nodes xp, Xi, x 2 m 
M. 2 . // rank [fi , f 2 , xi — xo,x 2 — xp] = 1 then there exists a purely elastic planar 
network with force response proportional to (f^, , fj) 2 ". T/ie internal nodes of 
such a network can be chosen to avoid a finite number of points and within an 
e— neighborhood of the convex hull of the terminals. 

Proof. We build up on the idea presented in 0, Theorem 5]. It can be easily 
observed (by choosing a coordinate system with Xp as origin) that for any point y 
in the plane, not collincar with xo, xi, x 2 , there is a force f , such that the following 
families, 

(21) ((f,y),(fp+f 2 -f,xp),(f 1 ,x 1 )) and ((-f , y), (f - f 2 , x ), (f 2 , x 2 )), 

form balanced systems of forces. Then, by using Lemma [7] there exists purely 
clastic networks with response matrices proportional with (f T , Lj + fj — f T , f^) T , 
and (— f T ,f T — fJjLj) 7, respectively. It is possible to choose the spring constants 
so that the associated quadratic forms are (see Remark [3]), 

q'(v, up, m) = (f • v + (f 2 - f) • u + f • up + L • ui) 2 and 

q"(v, up, u 2 ) = (-f • v + (f - f 2 ) • up + f 2 • u 2 ) 2 . 

Now, let us consider the infimum, 



(23) g r (u ,u 1 ,u 2 )= inf {q'(v,u ,Ui) + q"(v, u , u 2 )}. 

vSR 2 

The necessary condition of the infimum, 

(24) f • v = -1 ((f a - f) • up + f • u + L ■ Ul ) - \ ((f 2 - f ) • u - f 2 • u 2 ) , 
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implies the statement of the Lemma, i.e., 
(25) 



g(uo,Ui,u 2 ) = i (^>2 { * ' u ^ 



The point y and any additional points in Lemma [7] can be chosen to avoid a finite 
number of points in the plane, and within an e— neighborhood of the convex hull of 
the terminals. □ 

3.1.2. General rank one planar networks. We show in Lemma [TOl the construction 
of a network realizing a valid four terminal rank one response and then generalize 
the result to any number of terminals in Theorem [21 We start our argument with 
Lemma IS] which is a technical result needed later in this section. 

Lemma 9. Let fo,fi,f2,f3 and {xo, xi, x 2 , X3} be a balanced system of forces in 
K 2 . Then there exists a point y* in an e-neighborhood of the convex hull of the set 
{x ,xi,x 2 ,x 3 } such that y* ^ {x , Xi, x 2 , x 2 } and 

fi A Xi + fj A Xj - (ft + fj ■ ) A y* = 0, for some {i, j} e {0, 1, 2, 3}, i ^ j, 

or in other words, the three forces fi, fj and —(fi + fj) supported at the nodes Xj, 
Xj and y* are balanced. The point y* can be chosen to avoid a finite number of 
positions. 

Proof. For any pair of indices {i,j}, with i ^ j and i,j S {0, 1, 2, 3}, let ftj : K 2 — > R 
be defined by 

(26) fij (y) = fj A x, + fj A x,- - (ft + ft) A y . 
Using the balance of forces relations it can be easily observed that 

(27) fij = fji = -f kt = -f tk , for any {i,j,k,t} = {0, 1,2,3}. 
Let C be the convex hull of {xo, xi, x 2 , X3}, that is 

C = < z e R 2 J z = CjXj, for q > and Cj = 1 > , 

I i=0 i=0 J 

and C e be the e-neighborhood of the set C as defined in (TTTj) . 

Next we show that there exists a pair of indices {i,j} with i, j £ {0, 1,2,3}, with 
the property that there exists a point y* £ C e , such that ftj(y*) = 0. We reason 
by contradiction. Assume that the above is not true, i.e., 

For all pairs {i,j}, with i ^ j and i, j <E {0, 1, 2, 3} we have 

(28) f^ (y) ^0, for any ye C e . 

Using the continuity of the functions ftj and the convexity of the set C e , from 
(|28p we obtain that all the functions ftj have constant strictly positive or strictly 
negative sign over C £ . Using this observation, together with the relations (|2"7| . for 
any partition U {m, n} = {0, 1, 2, 3} we have 

(29) fy(y)f rml (y) < for y e C e . 

For simplicity we shall call from now the complement of ftj the function ft t with 
{k,t} = {0,l,2,3}\{i,j}. 

From ([2T| . (|28p and we conclude there are six different functions fy, with 
{i,i} G {0, 1, 2, 3}, out of which three functions have strictly positive sign on the 
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set C £ while their complements have strictly negative sign on C £ . Using this ob- 
servation, it can be easily checked that for {i,j,k,t} = {0, 1,2,3}, at least one of 
the following triplets of functions f^), (f i3 ,f lk ,{jk), (Ey, fit, fjt), (fij,fjt,fjk), 

(ffci, fjfe, fi t ), has a constant strict sign over the set C £ . Indeed, if f,j, f^ and fit do 
not have the same sign, cither two are positive and one is negative or two are nega- 
tive and one is positive. By replacing one by its complement we get three functions 
of the same sign. From the balance of forces relations one can immediately observe 
that the sum of the functions in any of these triplets is of the form 

(30) ± 2f l0 A (y - x, ) for some i £ {0, 1, 2, 3}. 

Finally relation (|30|) leads to a contradiction. Indeed for a triplet with the property 
that all the functions in the triplet have the same strict sign, the sum of its functions 
must have the same sign on C £ . However from (|30p the sum cannot have constant 
sign over the set C £ because it equals zero for y = x io . Thus the hypothesis ((28)) is 
false and we have that there exists a pair of distinct indices {i,j} C {0, 1,2,3} so 
that Ey(y*) = for some y* £ C £ . 

From ([30]) it is possible to choose y* ^ {xo,Xi,X2,X2}. Indeed from Lemma ITOl 
there are two indices {m, n} such that 

(31) fm«(y*)=o. 

Consider the point ys = y* + 5(f m + f n ) with 5 > 0. For S small enough it is clear 
that ys £ C £ \ {xo, Xi, X2, X3} and from (|2"faT) wc obtain 

f m n(ys) = fm A x m + f„ A x„ - (f m + f n ) A y s 

= fmn(y*) - S{f m + f n ) A (f m + f n ) 

which implies the statement of the Lemma. □ 

The next Lemma [10] shows that for any balanced system of four forces in R 2 
there exists a purely elastic four terminal planar network with proportional (rank 
one) force response. 

Lemma 10. Let fo, fi, f2, f3 and Xo, xi, X2, X3 be a balanced system of forces in M 2 . 
Then, there is a purely elastic four terminal planar network with force response 
proportional to (f^, f^, fj", f 3 r ) T . The internal nodes of such a network can be 
chosen away from a finite number of points, and within an e— neighborhood of the 
convex hull of the terminals. 

Proof. From LcmmalHlwe have that there exists a pair of indices {i, j}, a point y* £ 
C £ \ {x , xi, x 2 , x 3 }, and a force f = £j +fj, so that the sets (y*,f), (x,, fi), (xj, f,-) 
and (y„, — f), (xfe, ffc), (x t , ft) are balanced sets of forces. From the results of the 
previous section we have that both sets have a rank one network reproducing the 
forces. By rescaling the spring constants (see Remark [3]) their associated quadratic 
forms are, 

91 (v, Uj, Uj) = (f • v + fj • Uj + ij ■ uj) 2 
q 2 (v, u fc , u t ) = (-f ■ v + f fe • u fe + f t • u t ) 2 

where {k,t} = {0,1,2,3} \ {i,j}- The quadratic form for both networks taken 
together is 

g(uo,Ui,u 2 ,u 3 ) = inf{<7i(v, a,,^) + q 2 (v, u k , u t )}. 
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The optimality conditions are 



f • v = 



-(fj • U t + fj ■ Uj) + -(ffc ■ Ufc + f t • u t ) 



which yield the desired result 



) 



2 



q(uo, Ui, u 2 ,u 3 ) = - f 4 • ui 



□ 



The following Theorem is the main result of this section. We use the previous 
results for three and four terminal networks to prove the result in the general case 
of p-terminal networks by induction, following the approach of Camar-Eddine and 
Seppecher 

Theorem 2. Let f; and Xj, i = 1, . . . ,p be a balanced system of forces in R 2 . There 
is a purely elastic p terminal planar network with a force response proportional to 
(f^, . . . ,fp) T and with internal nodes in an e— neighborhood of the convex hull of 
the terminals and avoiding a finite number of points. 

Proof. We use an induction argument in the number of terminals as in 0, Theorem 
5], which we reproduce here for completeness. The p = 2 case corresponds to a 
single spring with same direction as the forces. The cases p = 3 and p = 4 are 
proved in Lemma [7J Lemma [8] and Lemma [TUl Assume for the induction argument 
that the theorem holds for any t < p terminals. Let r be the integer part of p/2, 
and let y be a node distinct from the terminals xq, ... ,x p . Then there are two 
forces f and f ' such that both families 



are balanced systems of forces, as can easily be seen by taking f = — (f + fi + f 2 + 
. . . + f r +i) and choosing Xi as the origin of coordinates. These families have r + 1 
and p — r + 1 terminal nodes, and both have less than p terminals when p > 4. 
By the induction hypothesis there are rank one networks with associated quadratic 
forms 



The quadratic form of both networks together is 

<f(ui,u 2 , . . . ,u p ) = infa'(v,ui, . . . ,u r+ i) + g"(v,u r+ i, . . . ,u p ). 

The optimality conditions are: 



(( f ,y), (fi + f',xi), (f2,x 2 ), . . . , (f r+ i,x r+ i)) and 
((-f,y), (-f',xi), (fr+i,x r+ i), . . . , (fp,X p )). 
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Thus the quadratic form of both networks is a rank one with 

1 / P 

g(ui, ...,u p ) = - I ^f, vu 

\i=i 

Notice that the additional point y can be picked inside an e— neighborhood of the 
convex hull of the terminal nodes and avoiding a finite number of points. Also we 
can use virtual crossings and trusses to make the network planar (see N2.3p . □ 

Remark 5. The networks in Lemmas^[RW[ \10l and Theorem^ may have crossing 
springs so are not strictly planar. To make them planar is suffices to convert all 
spring crossings with non-zero angle to a node as is done in [Q, Example 3.151. 
Zero-angle crossings can be eliminated by replacing springs with simple trusses ]Q, 
Example 3.2]. These network transformations are discussed in Wit. 31 

3.2. Networks in R 3 with rank one static response matrices. The construc- 
tion of these networks is essentially due to Camar-Eddine and Seppecher We 
first complete their construction of rank one four terminal networks to include some 
degenerate cases which correspond to planar networks (Lemma lll[) . Then following 
Camar-Eddine and Seppecher Q we use induction (Theorem [3|) to derive rank one 
networks with an arbitrary number of terminals. 

Lemma 11. Let fo, fi , f^, f 3 be a set of balanced forces at the nodes xo, xi, X2, 
X3 in R 3 . There is a purely elastic rank one network with force response propor- 
tional to ({q , , fJ) T . Moreover the internal nodes can be chosen within an 
e— neighborhood of the convex hull of the terminals and avoiding a finite number of 
points. 

Proof. As in the planar case there are two cases depending on the value of 

r = rank [f 1 ,f 2 ,f 3 ,x 1 -x ,x 2 -x ,x 3 -x ]. 

The construction for r = 3 is given in 0, Lemma 5]. When r < 2 the network is 
planar, so the result follows from Theorem [2] □ 

Theorem 3. Let fi and x^ ; i = l,...,p be a balanced system of forces in R 3 . 
There is a purely elastic p terminal network with a force response proportional to 
(ff, . . . , ip ) T . Moreover the internal nodes can be chosen within an e— neighborhood 
of the convex hull of the terminals and avoiding a finite number of points. 

Proof. The result follows from an induction argument similar to that of Camar- 
Eddine and Seppecher 0, Theorem 5] . See also the proof of Theorem [3] □ 

4. Characterization of the dynamic response function 

To fully characterize the dynamic response matrices, we take a function W(cj) 
as in Lemma H] and show that we can construct a network that has W(w) as its 
frequency response. The construction relies on the static case (Theorem Q]) and 
the existence of a network of springs and masses with rank one response that has 
exactly one prescribed resonance fLemma [T2|) . Both networks in R 3 and planar can 
be constructed. 

Theorem 4. Let W(cj) be a matrix valued function of to satisfying the properties 
of Lemma [^} Then for any choice of terminal node positions, there is a network 
( either planar or in M 3 ) of springs and masses with W(w) as its response function. 
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Moreover the internal nodes of such a network can be chosen to avoid a finite number 
of positions and within an e— neighborhood of the convex hull of the terminals. 

Proof. It is convenient to rewrite ijHJ) as, 

(32) W(w) = W(0) - w 2 M + V C«^-^ 2-. 

Since W(0) has the properties of Lemma [H by Theorem [1] there is a (static) net- 
work of springs that has W(0) as its response matrix. Since the are positive 
scmidefinite we can use the spectral decomposition to write 

rii 

C«=£A«cf(cff, for *=!,..., 

where the (Aj,Cj ) are the eigenpairs of CW with \y > and ni = rank(CW). 
By Lemma [T2] we can construct a network with response function 



, . (i) _i (») W /. (i) _i (i) 



T 2 

By the superposition principle there is a network with response the sum appearing 
in (f3"2"| . Finally to obtain the — w 2 M term, simply endow the terminal node Xj with 
a mass equal to the id— th diagonal element of M (The mass of node i is repeated 
d times in the nd x nd matrix M). To obtain a planar network, simply replace all 
spring crossings by "virtual crossings" (see ^2.3[) . □ 

4.1. Rank one response matrices with resonance. The following lemma shows 
how to design a network with arbitrary rank one response and one single resonance 
at a prescribed frequency. The network we construct has a purely dynamic response 
as it has a zero response matrix in the static case u> = 0. The construction is valid 
for both networks in Ml 3 and planar networks. 

Lemma 12. Let Xj and fj, i = 1, . . . , n, be arbitrary points and forces and ujq 
a given finite resonance frequency. There is a network with terminals X; composed 
of springs and masses with rank one response function 

W(w) = ff T ^%, where f T = (if , fj, . . . , fj). 

Moreover the internal nodes of such a network can be chosen to avoid a finite number 
of positions and within an e— neighborhood of the convex hull of the terminals. 

Proof. Let x„ + i and x„ + 2 be two distinct nodes and choose the forces f n +i and 
f n +2 so that the system (xj, f,), i = 1, . . . , n + 2 is balanced. Take for example a 
force f„ + 2 ^ in the line (d = 2) or plane (d = 3) 

(x n+2 - x„+i) A f„ +2 = - ^(Xj - x n+ i) A f;, 

i=l 

and choose f„+i such that J27=i ^ = 0. Then by Theorem [2] there is a rank one 
network with force response proportional to (f T , fj +1 , fJ +2 ) T - Attach a mass m 
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to nodes x„+i and ~x. n +2- The spring constants in the network can be rescaled 
(Remark [3]) so that the equations of motion are 

muj 2 \ij 

where a T = (^+1,^+2) 7^ Oj the displacements and uj are respectively at the 
"boundary" nodes xi,...,x„ and the "interior" nodes x n +i) x n+2- Then solving 
the system for the resulting forces wg at the "boundary" nodes we get, 

rp ( ||a|| 2 \ rp LO 2 

wb = ff T u B 1 + ^ 2 = ff T u B g . 

\ muj 2 — ||a|| y uj 2 — ||a|| /to 

Finally choose the mass m = ||a|| 2 /ujq. The position of the internal nodes x rl+ i 
and x n+ 2 is flexible and by Theorem [2] so is that of any interior nodes in the rank 
one network involved in the construction. □ 



Appendix A. Stability to small perturbations 

We show that the response function of an elastodynamic network is stable to 
changes in the network, which could come from cither modifying the spring con- 
stants of existing springs or possibly adding springs with small spring constants 
between any two nodes in the network. However we do not allow springs to be 
deleted from the network. We first show stability of the response of static networks 
and then stability for the response function of elastodynamic networks. 



A.l. Stability in the static case. Let A be the response matrix of an elastic 
network with all nodes considered as terminals and W be the response matrix at 
the terminal nodes as given by 0$. If we add or modify (but not delete) springs 
then the new response with all nodes considered as terminals is A + eE, e > 0, and 
its response at the terminals W(e). We prove the following result 

Lemma 13. Let e > 0. As e — > 0, we have W(e) — > W. 

The stability result for the static case may seem surprising at first because the 
pseudo-inverse we used to find the response at the terminals (U) is not continuous 
(see e.g. [1, §5.5.5]). However Lemma Q] guarantees that the instabilities are con- 
trolled as they remain (roughly speaking) in A/"(A//). Before showing Lemma [T3l we 
need to establish the following relation between the floppy modes of the perturbed 
and unperturbed stiffness matrices. 

Lemma 14. For all e > sufficiently small, Af(A]i + eE//) C Af(Aji). Moreover 
A/"(A// + eE//) is independent of e, and depends only on the connectivity of the new 
network. In other words if a network is perturbed by adding springs or modifying 
existing springs, then a floppy mode of the perturbed network must be a floppy mode 
of the unperturbed network. 
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Proof. If ui is a floppy mode of the perturbed network then 
= uj (A n + eE//)u/ = [0 uj] (A + cE) 



^2 ( k i,j + el ij) ( ( u * _ u 



old springs 



>i>.A", " . ! 



new springs 

iei. jeBui 



Since the perturbed network is a spring network, the new spring constants should 
be positive and all the terms in the sums above vanish, a condition which is indepen- 
dent of e. Since all the terms in the sums above vanish it follows that uj € JV(Ajj): 



Uj Aj/Uj = ^2 hj ( ( u * - u j) 



old springs 
i£l, jeBUl 



0. 



□ 



We are now ready to prove stability for the static case. 

Proof, (of Lemma [T3|) Let e > be sufficiently small. By Lemma [T4l it is possible to 
find a unitary matrix [U, V, W] independent of e such that 7?.([V, W]) = Af(Au) 
and 1Z(W) = N(Au + eE//). Writing An + e~Eu in the new basis gives, 



(33) 



A // + eE // = [U,V,W] 



A + eEi eE 2 







eE 3 




[U,V,W]^ 



where A = U T A//U. Because of our choice of basis both A and the non-zero 
block in (|33p must be invertible and symmetric positive definite. The inverse of 
this block is 



A + eEi 



eEj 



eE 2 
eE 3 



(34) 



A(e)- 1 (I + eE 2 E(e)- 1 E 3 - 1 E 2 n A(e)- 1 ) -A(e)- 1 E 2 E( e )- 1 E 3 - 1 



-E^E-^A^)- 1 



A" 1 + eGi + o(e) -A _1 E 2 E 3 _1 



-E^A" 1 



eG T 2 



o(e) 



l e; 



e- 1 E(e)" 1 E 3 

+ eG 2 + o(e) 
G 3 + o(l) 



where E(e) = (I - eE 3 " 1 E 2 r A(e)- 1 E 2 ) and A(e) = A + eEi. Notice that E 3 is 
a submatrix of a symmetric positive definite matrix and thus must be invertible. 
The second equality comes from the standard perturbation formula for the inverse 
(see e.g. [1, §2.3.4]) and the matrices Gi,G 2 ,G 3 are independent of e. Wc now 
examine the response of the perturbed matrix. Since the first term in ((4|) is linear 
in E, it is stable to perturbations. Now the negative of the second term in (j4|) can 
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be written as: 
(35) 

(A B i + eE B i)(A n + eE // ) t (A /B + eE IB ) = A B i(A n + eE H yA IB 
+ eE BI (A n + eEjj) t Aj B + eA BI (A n + eE n ^E IB + e 2 E BI (A u + eEjj^E/B. 
Moreover in the basis [U, V, W] the pseudo-inverse becomes, 

{Au + eEu^ = [U,V] 

Using Lemma[T]and LemmalMl we have A B j[V, W] = [0, 0] and (A B i+eE B i)W = 
0. The leading order asymptotics of each of the terms in (f35|) are 

A BI (A n + eE H yA IB = A B/ UA- 1 U T A /B + (9(e) = A BI A\jA IB + 0(e), 

eE BJ (A /7 + eE n )^A IB = eE B /(UA _1 U T - VE^ 1 E^A~ 1 U T )A /B + 0(e 2 ), 

eA B /(A// + cE/jjtE/B = eAs^UA- 1 ^ - UA- 1 E 2 E 3 - 1 V T )E /B + 0(e 2 ), 

e 2 E BI (A H + cEujtE/B = eE B/ VE 3 - 1 V T E J B + 0(e 2 ). 

which proves the desired result: 

(A BI + €E BI )(A n + eE // ) t (A /iJ + eE IB ) = A bi A\jA ib + <D(e). 

□ 

A. 2. Stability in the dynamic case. As in the static case, we deal only with 
network perturbations that modify existing springs or add new springs, but exclud- 
ing spring deletions. Denote by K the response when all the nodes are terminals 
and let e > be sufficiently small so that K + eE is the response of the perturbed 
network. We show the following result. 

Lemma 15. Partition the interior nodes I into massless nodes L and nodes with 

mass J, as in Lemma\^ Let u> be a frequency such that uj 2 is not an eigenvalue of 
r/2 — 1/2 

Mjj K,/,/M 7 j , and K + eE be the response of the perturbed network where we 
allow for new springs or changes in the spring constants, but no spring deletions. 
Then as e — ► 0, 

W(w;e) = W(w) + G(e), 
where W(cj;e) (resp. W(cu)) is the response at the terminal nodes of the perturbed 
(resp. unperturbed) network at frequency uj . 

Proof. By Lcmma[l3]tlie matrix K (the response matrix of the network with termi- 
nals BUJ, see ([7])) is stable to such spring perturbations, meaning that the response 
of the perturbed network at the nodes BUJ satisfies K(e) = K + eE + o(e), for some 
matrix E independent of e. Since both the perturbed and unperturbed responses 
are symmetric, the Wielandt-Hoffman theorem (see e.g. 0, §8.1.2]) implies that 
there is a reordering of the eigenvalues oj 2 (e) of K,/j(e) such that 

k 2 (e)- W f|<e||E|| F , 

with uj 2 being the eigenvalues of K,/j and || • || F denoting the Frobenius matrix 
norm. Therefore if e is sufficiently small and u> is not an eigenvalue of M 7 y KjjM j , 
then u is not an eigenvalue of Mjj^Kjj^Mjy 2 either and the matrix K(e) — 
w 2 Mjj is invertiblc. Thus using (|6|) and the standard perturbation formula for the 



A + eEi 

eE T 2 



eE 2 
eE 3 



[U,Vf 
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inverse, it is possible to show that the response at the terminals of the perturbed 
matrix is W(cj; e) = W(w) + 0(e). □ 

Appendix B. Eliminating floppy modes by adding springs 

We use the stability results from the previous Appendix to show that if a network 
has floppy modes then there is a network with no floppy modes and a response 
function arbitrarily close to that of the original network. Some examples of floppy 
modes (for planar networks) are given in Figure [3] Our strategy to remove floppy 
modes is to connect all nodes (be them terminal or interior nodes) of the network 
by springs with small spring constants, thus creating a complete graph with the 
nodes I U B and springs as edges. By Lemma [T3] the response of the new network 
can be made arbitrarily close to the response of the unperturbed network. 

Let A be the response of the network where all nodes arc terminals and let 
A + eE be response when we have added all these new springs. Then if uj is a 
floppy mode of the new network, proceeding as in Lemma [14] gives 



= uf (Ajj + eEu)uj = [0 uj] (A + eE) 



X ( k >.j + e-kj) l(ui -Uj) 



old springs 



x,- — X-, 



2 



ei<,j ( (Ui - ttr) 



v 3 



2 



new springs 
ieZ, jGBUI 

This is equivalent to saying that 

(36) Vt e I and j e BUI, (u ; - u.,) ■ (x t - Xj ) = 0. 

Assume that we are working in d dimensions and that we have d+l nodes Xi , . . . , x^+i 
that form a non-degenerate triangle (d = 2) or tetrahedron (d = 3) and for which 
Ui = U2 = ■ • ■ = Ud+i = 0. Then since every interior node y is connected to the 
xi, . . . , Xd+i, equation (|3"6")l implies that v • (y — x,) = 0, i = 1, . . . , d + 1, where 
v is the displacement associated with y. Since these special "anchor" nodes are 
non-degenerate, rank [xi — y, X2 — y, . . . , x^+i — y] = d and we must have v = 0. 
Repeating this for every interior node we get uj = 0, and so the network does not 
have any floppy modes. 

We now need to show which networks have "anchor" nodes. Clearly if the net- 
work has d+l terminal nodes forming a non-degenerate triangle (in d = 2) or 
tetrahedron (in d = 3), then the network does not have any floppy modes, since the 
terminal nodes do not move. 

If we are in d = 2 dimension and the terminal nodes do not form a non-degencrate 
triangle, then all terminals must lie on a line. Since the network has at least two 
terminal nodes (otherwise we cannot balance forces) , we can add two interior nodes 
with a truss as in H2. 31 without changing the response. Let xi,X2 be two terminal 
nodes and y be one of the interior nodes of the truss, with associated displacements 
Ui = U2 = and v. Then condition ([3^]) implies that v • (xi — y) = v • (X2 — y) = 0, 
i.e. v = 0. Thus the nodes xi,X2,y form an "anchor" and the network does not 
have any floppy modes. 
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Figure 3. Two types of floppy modes. The terminals are white 
circles and the interior nodes are in black. The direction in which 
the node can move with zero force is given with arrows. 



If we are in d — 3 dimension and we cannot form a non-degenerate tetrahedron 
from the terminal nodes, then the terminals must lie on a plane. Assume further 
that the terminals do not lie on a line, we shall deal with this case later. Let 
Xi,X2,X3 be three terminal nodes forming a triangle. Then replacing e.g. the 
spring between Xi and X2 by a truss (as in ij2.3[) . we introduce three new interior 
nodes and at least one of them y is not in the plane where Xi,X2,X3 lie. Let v 
be the displacement of y and Ui = U2 = 113 = be the displacements of the 
boundary nodes. Then condition (|36|) implies that v • (xj — y) = for i = 1,2,3. 
Since rank [xi — y, X2 — y, X3 — y] = 3, we must have v = 0. Therefore the nodes 
Xi,X2,X3,y form an "anchor" and the network does not have any floppy modes. 

If we are in d = 3 dimension and all the terminal nodes lie on a line then 
every interior node y forms a triangle with two terminal nodes Xi , X2 . In this case 
condition |36|) means that the displacement v of node y is orthogonal to the plane 
formed by xi,X2,y, and in particular to the axis where all terminals lie. Thus 
in this case the floppy modes cannot be eliminated by adding springs or interior 
nodes, as any additional interior node is in this situation as well. This corresponds 
to rotations of the network around the axis where all the terminals lie. 
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